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Abstract 

In the present work, we investigate the sohton structure of the Vakhnenko equation with a 
'dissipative' term, by means of the Hirota's method. As a result, we unearth three kinds of sohton 
solutions depending upon the dissipation parameter. We further find that the scattering behavior 
among such structures exhibits a peculiar feature that may be the fission phenomenon. 
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I. INTRODUCTION 



Nonlinear systems and their related waves in different branches of physics have recently 
been investigated in detail through theory and experiments. Since the investigation of the 
nonlinear transmission line, pioneered by Hirota and Sazuki various, kinds of nonlinear 

rinnn 

wave propagation have been discussed [2|, y, 14, 15[. Integrable models of such systems have 
3een applied to many fields of research such as condensed matter , fluid mechanics 
8L plasma physics (9|, optics 10, 11|, communications 3], chemistry [l^, and biology 



151 . Il6l | . These integrable models which may take the form of nonlinear partial differential 
equations, may possess some kinds of self-confined solitary waves solutions known as solitons. 

Since Gardner, Greene, Kruskal and Miura have solved the Korteweg-de Vries equation 
using the inverse scattering method [l7|, the modern theory of solitons has been deeply 
studied in mathematics, and widely applied in physics. A wealth of interesting properties of 
soliton equations have been found . Indeed, a soliton equation may possess Backlund and 
Darboux transformations, the Lax-pair, N-soliton solutions, a bilinear form, a zero-curvature 
form, the Painleve property, infinitely many symmetries and conservation laws, just to name 



a few 



18|. 



Recently, Vakhnenko [19| have investigated high-frequency perturbations in a relaxing 
barothropic medium. He has then derived a new nonlinear evolution equation known as 
Vakhnenko equation |l^ given by 

utx + {uux)^ + aux + M = 0, (1) 

where u may stand for a physical observable, and subscripts t and x appended to u may 
refer to partial differentiation with respect to these independent variables. It may be noted 



here that the aMa;-term may be regarded and defined as the 'dissipation' term 
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JJL Without 
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'dissipation', this equation which has been subject to many investigations 
may possess loop-like soliton solutions. As far as we are concerned, few investigations 
of Eq. ([1]) may Kave been done. Va.Knen.o Q Kas pointed on. E, (LLJ) may possess 
ambiguous solutions in the form of a solitary wave, and he has proved that the dissipative 
term, with a dissipative parameter a less than that limit value, does not destroy the loop-like 
solutions. One underlying query that may obviously be asked may be related to the existence 
of other kinds of solutions when the dissipation parameter is greater or equal to this limit 
value. Besides, the question of soliton structure of such solutions may be investigated. 



In order to provide some answers to these queries, the study of the fundamental role 
played by the dissipation parameter may be done around the following points. In Sec. [Tll 
the different kinds of the one-soliton solution are presented. In Sec. Illlt the two-soliton 
solutions and their scattering behavior are depicted. Some peculiar phenomenon such as the 
fission phenomenon, has been identified. Sec. IIVI is devoted to a brief summary. 

II. BILINEARIZATION OF THE VAKHNENKO EQUATION AND CONSTRUC- 
TION OF THE ONE-SOLITON SOLUTIONS 



Introducing new independent variables X and T as follows [22 1 



x = T+ U{s,T)ds + xo, t = X, (2) 

J — oo 

where u(t,x) = U (X, T) , and xq stands for an arbitrary constant, setting U = Wx, Eq. ([1]) 



transformed to 



22| 



WxxT + WxWt + aWr + Wx = 0. (3) 

By taking 

W = 6i\nF)^, (4) 
Eq. ([3]) may take the following coupled bilinearized forms 

{DtD^x + D\) F-F + aGF = 0, (5a) 

Dt {D\F -F) -F^ - GF^ = 0, (5b) 

where G may stand for an arbitrary function. The quantities Dt and Dx denote the Hirota 
operators. According to the usual procedure, the soliton solutions may be constructed by 
expanding F and G in suitable formal power series. 
Thus, the one-soliton solutions may be given by 

F = l + exp{2ri), (6) 

where 

7] = KX-u;T + r]o, (7) 
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quantities K and uj standing for wave number and angular frequency, respectively. The 
dispersion equation may be given by 

2 (a + 2ir) - 1 = 0, (8) 

and uj = Kv, v > being the velocity of the wave. Combining Eqs. (jlj) and may lead to 

W = 6K[1 + tanh (r^)] , (9) 

and hence 

U = f/Msech^ (ri) , (10) 

where the amplitude Um = QK"^- 

In order to discuss the different types of solutions, it seems worthy to consider the fol- 
lowing equation 

dT=[l- QKusech^ (r/)] d^. (11) 

Thus, setting A = QKuj, we may distinguish three kinds of solutions according to the cases 
A < 1, A = 1, and A > 1. As a result, we find that 

1. for a < loop-like solution may be obtained (see FIG. 1 where loop is represented 
by the solid line); 

2. for a = cusp-like solution may be obtained (see FIG. 1 where cusp is represented 
by the broken line); 

3. finally, for a > hump-like solution may be obtained (see FIG. 1 where hump is 
represented by the dotted line). 

In FIG. 1, the previous kinds of soliton solutions have been depicted as variations of U /Um 
vs X. It seems also worth depicting the domain in which these solutions may exist. In this 
view, we depict at FIG. 2, the variations of a vs the velocity of waves v. In this figure, the 
white area may refer to hump-like soliton solutions. It seems also worthy to compare the 
amplitudes of these solutions. The variations of Um vs the velocity of waves v are depicted 
in FIG. 3 where the white area may refer to the loop-like soliton solutions. Thus, for a 
given velocity, the greatest amplitude may be due to loop-like solutions and the smallest 
amplitude belongs to hump-like solutions. 



III. THE TWO-SOLITON SOLUTIONS 

The solution to Eq. ([5]) corresponding to two-soliton solutions may be given by 

F = 1 + exp(r/2) + exp(2?72) + Aexp{Arii) + B exp{Ar]2) + Cexp2{r]i + r/2), (12) 
where the phases rji {i = 1,2) may be given by 

rji = KiX - UiT + rjio, (i = l,2), (13) 

and the dispersion relations 

2{a + 2Ki)ui-l = 0, {i = l,2) (14) 

where iVi = KiVi {i = 1,2), vt standing for velocities of waves. 
The coefficients A, B and C may be given by 

2[a + bKi) 
2 (« + 6A2) 

^_ 2a [(cui - UJ2) {Kl - Kl) + 2K^K2 {uji + UJ2)] + 4 (cu2 - cuQ {K^ - K^f + {K^ - K^f 

{K^ + K^f [2 {uJr + UJ2) {a + 2K^ + K^) - 1] 

(15c) 

In order to investigate the scattering behavior among these two-soliton solutions, we may 
firstly consider some underlying snapshots. 

• For a = 1.2, according to Eq. (fTTl) . two snapshots depicting loop-like and hump-like 
solitons may be depicted at FIGS. 4 and 5 at times t = —15 and t = 11, respectively; 

• for a = 0.1, the same phenomenon may be observed but we may see that the loop-like 
soliton splits into two other loops of the same amplitudes. This case may be depicted 
at FIGS. 6 and 7 at times t = — 15 and t = 11, respectively. Thus, as a result, there 
may exist a value at which the loop-like soliton begins to split. It may be suggested 
that for velocities Vi = 0.24 and V2 = 0.12, there may not be splitting provided the 
relation 1.2 < a < 2.6 holds, whereas for a < 1.2, splitting may be observed; 



• for a = 2.6 with velocities vi = 0.24 and V2 = 0.12, according to Eq. (ITTi) . two 
snapshots depicting cusp-fike and hump-fike solitons may be depicted at FIGS. 8 and 
9 at times t = — 15 and t = 11, respectively; 



• finally for a = 5.0 with velocities vi = 0.24 and V2 = 0.12, according to Eq. (fTTI) . two 
snapshots depicting large hump-like and small hump-like solitons may be depicted at 
FIGS. 10 and 11 at times t = —15 and t = 11, respectively. 

Now, the full scattering behavior among these soliton solutions may be properly analyzed. As 
it may be observed, for negative values of time t, only one-soliton solution may be observed. 
For positive values of time t, the initial one-soliton solution splits into two other solitons 
with an increase in the amplitude. The soliton with the small amplitude may have hump- like 
shape. This kind of phenomenon that may refer to a fission, has recently been investigated 
by Morrison and Parkes while deriving the N-soliton solution to a generalized Vakhnenko 
equation. 



IV. SUMMARY 



In this paper, we have investigated the soliton structure of the Vakhnenko equation 22 



with a 'dissipative' term, by means of the Hirota's method. This model equation derived 
by Vakhnenko [2^ may be underlying in description of high-frequency perturbations in a 
relaxing barothropic medium. As a result, we have unearthed three kinds of soliton solutions 
depending upon the dissipation parameter. These soliton solutions may have loop-like, cusp- 
like and hump-like shapes, respectively. We have further found that the scattering behavior 
among such structures exhibits a peculiar feature known as fission. Indeed, the initial one- 
soliton splits into two other kinds of solitons as time elapses from negative values to positive 
values. This kind of phenomenon may have been pointed out by Morrison and Parkes 25 1 
while investigating the N-soliton solution to a generalized Vakhnenko equation. In order 
to provide the full detail on the interactions such as asymptotic bahavior of the scattering 



among the soliton so 
described in Ref. 



22| (see Eq. ([T])), the same procedure 



utions to the Vakhnenko equation 
25l ] may be followed. Nonetheless, for some convenience due to the length 
of the paper, we do not report such result here. We may have only focused our interest to 
the effect of the dissipative parameter a on the Vakhnenko equation [22i] (see Eq. ([T])). 
For a further interest, it may be worthy to investigate the N-soliton solutions (A^ > 3) 



to the Vakhnenko equation 22|] (see Eq. ([T])). Another interesting study may be to fix the 
dissipation parameter a, and subsequently investigate the effect of the velocities of waves on 
the system. Besides, following a recent work of Konno and Kakuhata ^] on rotating soliton 



solutions to a coupled dispersionless system, it may be worth considering complex-valued 
angular frequency uj and wave number K into the dispersion equation. This may hopefully 
help investigating rotating soliton solutions to the Vakhnenko equation |22j] (see Eq. ([T])). 
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FIG. 2: Variations of the dissipation parameter a vs the velocity v of waves. White zone may refer 
to hump-hke soUton solutions. 
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FIG. 3: Variations of the amplitude Um of waves vs the velocity v of waves. White area may refer 
to loop-like soliton solutions. 
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FIG. 4: Snapshot of loop-like soliton solution at time t = —15 from Eq. (jl2p for a = 1.2, vi = 0.24 
and V2 = 0.12. 
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FIG. 5: Snapshot of loop-like soliton solution at time t = 11 from Eq. ()12p for a = 1.2, vi = 0.24 
and V2 = 0.12. 
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FIG. 6: Snapshot of splitting loop-like soliton solution at time t = —15 from Eq. (jl2p for a = 0.1, 
VI = 0.24 and V2 = 0.12. 
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FIG. 7: Snapshot of splitting loop- like soliton solution at time t = 11 from Eq. (jl2p for 
VI = 0.24 and V2 = 0.12. 
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FIG. 8: Snapshot of cusp-like soliton solution at time t = —15 from Eq. (jl2p for a = 2.6, vi = 0.24 
and V2 = 0.12. 



16 




FIG. 9: Snapshot of cusp-like soliton solution at time t = 11 from Eq. (jl2p for a = 2.6, vi = 0.24 
and V2 = 0.12. 
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FIG. 10: Snapshot of hump-like sohton solution at time t = —15 from Eq. (|12p for 
VI = 0.24 and V2 = 0.12. 



18 





-800 -600 -400 -200 200 400 



FIG. 11: Snapshot of hump-like soliton solution at time t = 11 from Eq. (jl2p for a = 5.0, vi = 0.24 
and V2 = 0.12. 
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